A differential equation approach for examining the subtraction schemes 
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We propose a natural differential equation with respect to mass(es) to analyze the scheme de- 
pendence problem. It is shown that the vertex functions subtracted at an arbitrary Euclidean 
momentum (MOM) do not satisfy such differential equations, as extra unphysical mass dependence 
is introduced which is shown to lead to the violation of the canonical form of the Slavnov- Taylor 
identities, a notorious fact with MOM schemes. By the way, the traditional advantage of MOM 
schemes in decoupling issue is shown to be lost in the context of Callan-Symanzik equations. 
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Conventionally, it is widely believed that different 
renormalization schemes (that differ arbitrarily with each 
other) are perfectly equivalent provided the perturbation 
is completed. Once truncated at finite order, they yield 
the annoying scheme dependence problem i.e., differ- 
ent schemes yield different predictions, and we have no 
better reasons to prefer one scheme to the others. But 
recently, it has been realized that the traditional on-shell 
definition of unstable particles masses and width (say, 
e.g., masses and widths of W^,Z^ and Higgs particles 
in the electroweak model) is problematic and should be 
replaced by the more physical and appropriate pole mass 
schemes [|], as the on-shell schemes would incur severe 
gauge dependence and IR singularities or ambiguities. 
That means, though we do not know the unique way of 
performing the physical definition, the freedom in choos- 
ing renormalization conditions or schemes deserves closer 
examination. 

In this brief report, we suggest to examine the prob- 
lem by formulating the Feynman amplitudes as a series 
of differential equations parametrized in terms of mo- 
menta (and masses for certain theories)^, this is feasible 
since the loop amplitudes must be functions of exter- 
nal momenta and masses after all. Then, as we shall 
see shortly, not all subtraction prescriptions or schemes 
satisfy such differential equations. The differential equa- 
tion analysis follows from the well known fact that the 
differentiation with respect to physical parameters like 
momenta and masses lowers the divergence degree of a 
superficially divergent Feynman diagram Q]. The solu- 
tions to these natural differential equations must be well 
defined or finite functions of momenta (and masses if any) 
and uniquely defined up to a certain set of constants to be 
defined through 'boundary' conditions, which just corre- 



*Viewing the quantum field theories as just effective sectors 
of the complete and well defined quantum theory of every- 
thing (QTOE), we can derive finite loop integrals by formu- 
lating them as solutions of appropriate differential equations 
in terms of energy, momenta and masses, without introducing 
any UV divergences, see Ref. 



spond to the definition of a renormalization/subtraction 
scheme, the parameters of the solutions must have been 
finite ones. It is worthwhile to point out that such analy- 
sis has already applied in the usual renormalization pro- 
grams through the critical use of Ward identities, a set 
of partial differentiation equations in terms of momenta. 
Our differential equation analysis is in fact a mathemati- 
cal generalization of the differential form Ward identities. 
In this short report, we temporarily focus on the differen- 
tial equations in terms of masses for QFTs with massive 
fields. 

One might object such analysis with respect to masses 
as the differentiation operation would shift or perturb the 
original theory. Such an objection is simply unjustified. 
First, no matter what kind of changes the differentia- 
tion brings into the theory, it does bring us information 
between the two theories that slightly differ in masses. 
Then as long as the difference is essentially analytical, we 
can integrate the difference back to find mass dependence 
of the original theory (up to ambiguities to be fixed by 
boundary conditions). Second, a more direct and math- 
ematical way to see this is, if we believe that the theories 
depend on the mass parameters in an almost analytical 
way, or that the Feynman diagrams are regular functions 
of masses, then no matter what the differentiation with 
respect to masses physically brings about, we can per- 
form differential analysis. Until now, we have not learned 
of any example of QFT whose mass dependence is so sin- 
gular that it defies differentiation equation analysis. So, 
as far as the mass dependence of Feynman amplitudes is 
concerned, it does not matter whether or not a method 
shifted the original theory or violated any symmetry of 
the original theory as long as the perturbation or viola- 
tion give rise to no completely singular mass dependence. 
We do not know any special roles played by the gauge 
symmetry and SUSY in preserving the analytical mass 
dependence. 

In fact such operation preserves most novel symme- 
tries, like gauge invariance, Lorentz invariance, and 
SUSY since canonical masses are gauge invariant, 
Lorentz invariant and SUSY invariant if they are masses 
of SUSY multiplet. For masses from symmetry breaking. 
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the difFerentiation equation analysis will become compli- 
cated due to the complexity of the Higgs sector which will 
be examined in a separate report. Here we focus on the 
invariant mass cases to demonstrate that not all subtrac- 
tion schemes are consistent with the differential equation 
analysis. Note that we are not claiming that these sub- 
traction schemes can not consistently renormalize QFTs 
with masses. We leave the implications of our result open 
to the readers. 

We will illustrate the analysis with a simple vertex 
function at the lowest order in a QFT, say, the 1- 
loop photon vacuum polarization tensor H'^'^ {p, — p, m) 
(= -ie^ J d"ktr{-f^'JJ^-t''J^}) in QED for simplic- 
ity. But the main points of the analysis apply in the same 
way to higher orders and to other theories like massive 
scalar theories and heavy quark theories. 

It is easy to see that this amplitude satisfies the fol- 
lowing well defined inhomogeneous differential equation 
in any gauge invariant regularization scheme (GIR) 
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as the RHS of Eq.(l) is well defined in such schemes. 
Or it satisfies the following well defined equation in any 
regularization schemes 
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where Cf is the combinatorial factor arising from the 
difFerentiation operation. 

It suffices to demonstrate our points with Eq.(l). 
Dropping the mass independent transverse factor 
{g^'^p^ — p^p'^) we have the following equation 



d™n(p , m) = — / dx— — 

-^^ /q m-^ — x{l — X)p^ 



(3) 



The solution to this equation is easy to find, it reads 



m? — x(l — x)p'^ 
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with C being the natural integration constant indepen- 
dent of to be fixed by regularization and/or subtraction 
schemes and F (p^) denoting an arbitrary mass indepen- 
dent function of momentum which can be fixed to be zero 
by solving the differential equation in terms of momen- 
tum. Then we have 
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with mass here being a finite or renormalized one. Note 
that the constant C here denotes any regularization 
scheme that is both gauge invariant and mass indepen- 
dent. In other words, those possible mass dependent reg- 
ularization schemes are already excluded, just like we ex- 
clude those regularization schemes that violate gauge in- 
variance and other physical symmetries in realistic mod- 
els. 

In all mass independent schemes, the subtraction is so 
defined that C is fixed to be a mass independent con- 
stant which still satisfies Eq.(3). However, in the MOM 
schemes, the subtraction is defined in the following way. 
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which ceases to be a solution of Eq.(3) but satisfies the 
following equation instead 
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It is clear in 



with 5(n^,m) — dx-^. — T-^ ^ 

Eq.(7) that the canonical relation of Eq.(3) is violated 
by the presence of an extra term that is dependent on the 
arbitrary subtraction point pi when massive particles are 
present. The extra term comes from the nontrivial depen- 
dence of subtraction constant upon particle mass in con- 
trast to the mass independent schemes (MIS) where 
dmCMis = 0, and therefore becomes the source of prob- 
lems in realistic applications. Regarding the mass depen- 
dence, the on-shell scheme belongs to the same class as 
MOM. It will lead to 'extra unphysicaP mass dependence 
for the vertex functions (including self-energies) for theo- 
ries with massive bosons, which is just the source of trou- 
bles illustrated in The MOM like schemes are also 
known to violate the canonical Ward or Slavnov- Taylor 
identities except in the background gauge |^], which we 
think is related to the fact pointed out above. 

We would like to provide more arguments regarding 
the mass dependence of renormalization schemes. Since 
the mass differentiation would insert mass operators that 
is closely related to the trace of the energy stress tensor 
which in turn couples to gravity, such operation should 
not lead to any 'extra unphysical' piece. It is easy to see 
that in Eq.(3) we can replace dm by mdm without chang- 
ing the problem. Alternatively, the mdm operation is 
part of the overall rescaling operation which also includes 
p^dpi-i . The latter operation should unveil the physi- 
cal dependence of vertex functions on spacetime, and 
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hence the mass dependence of vertex functions should 
also be physical due to the overall rescaling entangle the 
mass dependence with that of momenta. We can then 
understand why mass dependent schemes often violate 
the canonical Ward identities, as the differential form of 
Ward identities contain the operation dpi^, if we multi- 
ply this operator from the left by then we see that it 
is part of the scaling equations. Since the overall homo- 
geneity of the vertex functions with respect to all massive 
variables will entangle the momentum dependence with 
the mass dependence, the extra term in Eq.(7) will lead 
to violation of the canonical Ward identities. In this con- 
nection we note that the pole mass for fermion is related 
to the effective mass which is gauge invariant and IR fi- 
nite by the equation me//(p^, m)|p2^^,/2 = Mpoi 
[H, i.e., the dependence of pole mass on the MS running 
mass (Lagrangian mass) is normal since me/ / can be de- 
fined in MS scheme though it is scheme independent 
The situation for bosons is similar, see 

Thus, contrary to the naive expectation, the differen- 
tial equation analysis does not allow for all kinds of sub- 
traction schemes. That means we should reexamine our 
conventional belief with due care. If we accepted this ap- 
proach, then the conventionally claimed scheme invari- 
ance of certain objects within a subset of schemes 
(i.e., within the mass independent schemes) can now be 
viewed as a full 'invariance' for all schemes that solve 
the differential equations in terms of all physical param- 
eters. Otherwise, the significance of the claims like the 
significance of scheme invariance of the first two coeffi- 
cients (/3o,/3i) of beta function would be diminished by 
the mass dependent schemes. In a sense we advanced a 
plausible rationale for using mass independent schemes in 
our calculations. We also wish to point out that our in- 
vestigation in fact derives rationale from the well known 
fact that not all regularization schemes would satisfy 
the Ward identities. In our point of view, complicated 
schemes (like MOM) not only bring about inconvenience, 
but also might be unnatural choices. In other words, the 
freedom in choosing the schemes is a restricted freedom if 
one works in the differential equation approach suggested 
above. 

Indeed there exists no literature up to date that pro- 
vided convincing constraints on mass dependence of 
schemes, but it is also hard to convince people that a 
scheme that defines the one loop vacuum polarization 
tensor like below would be a physically acceptable choice: 

n^''(p, -P.-m; ^j)\\Toy 
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that is to say, such a subtraction leaving a finite term like 



the belief that the subtraction can be arbitrary as long 
as divergences are removed, at least the anomalous di- 
mension will be awkward-looking and very different from 
the standard results. In fact we can consider the prob- 
lem from the reverse angle, i.e., suppose we found the true 
physical and hence scheme independent parametrization, 
these physical parameters (together with a RG and renor- 
malization scheme invariant scale h-phys ) work coherently 
to give the unique mathematical formulas that are often 
complicated functions. Then it is impossible to trans- 
form these physical parameters into arbitrary things be- 
yond novel symmetries or invariances, i.e., it is impossi- 
ble to redefine the parameters in arbitrary schemes, since 
we must confront our calculations with experiments and 
experiments will only justify the schemes in which the 
dependence of the measured quantities (defined as func- 
tionals of Feynman amplitudes) upon the renormalized 
parameters are not quite different from the dependence 
upon physical parameters, as the renormalized parame- 
ters functionally take the same role as the physical ones, 
an obvious fact in the formulation of QFT. Our inves- 
tigation above seems to add doubts to the conventional 
attitude to this issue. Since it is just a first attempt, we 
will refrain from claim anything, we only wish to draw 
attention to the above aspects. 

The differential equations can be generalized to any 
order and any QFT with massive fields: 



(5™^)"r(pi,-..,(TOf);[q) 
r,. (pi,[0L,-..,(r<);[c]'), 



1,2, u>l, (9) 



( u^i+m^ )^°°° '^^^ '^'^^ acceptable choice, contrary to 



where (mf) refers to the various masses in the theory 
with s = 1 for fermion and s = 2 for boson, [C] and 
[C]' refer to the constants that will appear in the solu- 
tions and Lo being any times that will annihilate some 
of the constants, which is the secret of the analysis in 
mass differentiation. Note that we have made use of the 
fact that the differentiation with respect to mass(es) in- 
serts the mass vertex operators-the operator ((/j^/s) (for 
fermion it should be understood as 'ip'4')- 

Conventionally the MOM schemes are held to be ad- 
vantageous over the mass independent schemes in ex- 
hibiting good decoupling behavior required by physics. 
In this connection, we wish to prove that in the con- 
text of Callan-Symanzik equation |l|] the decoupling of 
heavy fields is achieved in the same way in both 
mass independent schemes and mass dependent schemes. 
Especially, we do not need the so-called "effective field 
theories" framework [|l^ to help the mass independent 
schemes, which is needed in the context of rcnormaliza- 
tion group equation (RGE) since it does not account for 
the effects of all the dimensional parameters. Such a 
proof is not seen yet in literature. Again we will illustrate 
it with a simple model, QED with a massive fermion in 
addition to ni massless fermions. In a mass independent 
scheme the Callan-Symanzik equation reads. 
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{Xdx - f3ada + 7r - Dr}T{{Xp),m,a, ^i) 
-«r®((Ap),TO,a,/i) 



(10) 



where = [1 + jmj'mip'ip, /3, 7r and are mass in- 
dependent functions of the renormaHzed coupHng a and 
all quantities are renormalized ones. At lowest order, 

When the mass goes to infinity, it is natural to expect 
that 

- zr®((p),TO, a,^)|„_oo 

= {A(3ada-Ajr)T{{p),a,fi), (11) 
-> {Xdx - Puada + iT-u - Dr}r((Ap), a, /i)|/.s. = 0. 

(12) 

Here (iis = (i + Af3, jnis = 7r + Ajr with the delta 
contributions coming from the mass insertion part in the 
infinite mass limit that will cancel the heavy fields' pieces 
in P and 7. The generalization to other theories with bo- 
son masses is an easy exercise. From Eq.(12) we see that 
the decoupling of heavy particles is realized in a natural 
way in the contexts of Callan-Symanzik equations. 

To verify the above deduction it is enough to demon- 
strate Eq.(ll) at the lowest order which is closely related 
to the observation that heavy particle limit provides a 
convenient algorithm for calculating trace anomalies fl^ ] 

2a 
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m(l -|-7m)V'V'|| 



(13) 



with EE -ie^-ff'ip and = = A-fA- When 

translated into Callan-Symanzik equations Eq.(13) is just 
Eq.(ll). The cancellation of the heavy particle contribu- 
tions is obvious since /3 -|- A/3 = + 1] — |^ = |^ri;. 

While in the MOM schemes the Callan-Symanzik equa- 
tion reads. 



{ Xd\ - fiMomCtda + 7r:Afo 

= -«rMo™((Ap),TO,a,/i) 



-Dr}rAfo,„((Ap), TO, a, ^l) 
(14) 



with the beta function, etc. being defined as Pmoiu = 
{fidf_i +m(l + )d,n)a, • • • , in contrast to the RGE 

definition: /J^/^^ = iJ-df^a, • • • , due to the mass depen- 
dence of the renormalization constants^. The resulting 
PMom also exhibits nondecoupling feature as in the mass 



^In the standard form of Callan-Symanzik equation, the 
mass operators inserted vertex functions appear separately 
and the definitions of (3, 7 in the MOM like schemes must in- 
clude m(l -I- ')m;Mom)dm to accouut for the 'full running' of 
the relevant renormalization constants. But in the alternative 
form of Callan-Symanzik equation (CSE), i.e., {Xdx — Pada + 
m{l + 7m)9m + 7r — -Dr}r((Ap), m, a, n) — 0, the definitions 



independent schemes. For example, at the lowest order, 
from the definition given above, a heavy particle's con- 
tribution to /3 at lowest order is 



/3. 
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Then it is easy to find what is similar to Eq.(ll) 
{{p),m,a, iJ,)\ 

= {AfiMomOlda - A'yr-Mom)'^Mo7n{{p),a.,y)\ls, (16) 

with APMom = -|f as Eq.(13) is also true in the MOM 
schemes at the lowest order. It is known that the first 
loop order beta of RGE in mass independent schemes 
differs from that in the MOM schemes Q. While in the 
context of Callan-Symanzik equation the first loop order 
beta function is the same in all schemes at one loop order 
(C.f. Eq.(15)), which in turn implies that the same de- 
coupling mechanism works in MOM schemes in the con- 
text of Callan-Symanzik equation. Thus the advantage of 
MOM schemes over the mass independent schemes in the 
decoupling issue is lost in the context of Callan-Symanzik 
equation and the decoupling is realized in the same way 
specified by Eq.(ll). 

We should note that, many MOM scheme calculations 
of couplings ijl^ are in fact calculations of the effective 
couplings with the momentum transfer defined at an Eu- 
clidean point. The CWZ [|l5| scheme for calculations in 
the presence of heavy quarks is in fact a MS scheme with 
the running scale fixed to be the particle mass, since in 
this scheme the mass also serve as a running scale for 
massive loops and the /3 function, etc. are all mass in- 
dependent. We should also note that our analysis only 
shows that the mass dependent schemes like MOM and 
on-shell schemes might cause problems for massive and 
some related sectors (like self-energy vertices, etc.) in a 
QFT. For massless sectors or massless QFTs, the prob- 
lem does not materialize. 

If one adopts our proposal for renormalization men- 
tioned above 0] , then the decoupling of heavy particles 



of the l3, 7 etc. are the same as in RGE. In this connection, we 
note that unfortunately many people have confused CSE with 
RGE which are different things in principle, CSE describes full 
scaling behavior while RGE only states finite renormalization 
effects. In addition, in the standard form of CSE, the trace 
of stress tensor is singled out, which describes the physical 
coupling of matters to gravity and therefore should be inde- 
pendent of renormalization schemes. Then we can anticipate 
that the heavy particle limit of this term should leave the 
same efi'ects to light sectors in any renormalization scheme, 
i.e., the same cancelling contributions to the full /3, 7 etc., just 
as what we will demonstrate here. 
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can be achieved in an automatic way as our proposal 
stands on the philosophical ground of the effective field 
theories please see Ref. Jf6[ . 

In summary, we suggested a differential equation anal- 
ysis of the radiative corrections in terms of momenta 
and masses and demonstrated that only certain subtrac- 
tion schemes satisfy such equations. Some conventional 
schemes like MOM failed to solve these equations. By the 
way, we gave a demonstration that in the context of the 
Callan-Symanzik equations all the schemes facilitate the 
decoupling of heavy particles in the same way, weakening 
the argument that the MOM like schemes is superior to 
mass independent schemes in decoupling issue. 
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